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We study theoretically the low-energy electric-dipole response of Y" = 20 isotones. We present 
results from a quasiparticle random-phase approximation (QRPA) and a continuum random-phase 
approximation (CRPA), and we compare them with results for the mirror Z — 20 nuclei. According 
to our analysis, enhanced El strength is expected energetically well below the giant dipole resonance 
in the proton-rich isotones. Large amounts of El strength in the asymmetric Y = 20 isotones are 
predicted, unlike their equally asymmetric Z — 20 mirror nuclei, pointing unambiguously to the 
role of structural effects such as loose binding. A proton-skin oscillation could develop especially 
in "^^Fe. The proper description of non localized threshold transitions and the nucleon effective 
mass in mean-field treatments may affect theoretical predictions. We call for systematic theoretical 
investigations to quantify the role bulk-matter properties, in anticipation of measurements of El 
transitions in proton-rich nuclei. 


I. INTRODUCTION 

The low-energy electric-dipole transitions of stable and 
unstable nuclei have been under investigation for more 
than two decades. The effort is posed to continue, as 
new experiments on exotic nuclei are being analyzed and 
new facilities with improved capabilities are being devel¬ 
oped. In the future the theoretical feedback may include 
results from methods drawing on ab initio nuclear struc¬ 
ture, namely the use of microscopic nuclear interactions. 
The concerted effort promises not only valuable empir¬ 
ical input into nucleosynthesis studies, but also crucial 
discoveries regarding nuclear structure far from the sta¬ 
bility valley, and the nuclear equation of state. 

The origin of low-energy El strength remains under 
debate, even in the case of stable nuclei. A particularly 
vexing question is whether or not highly asymmetric nu¬ 
clei, with an excess of neutrons (or protons), can develop 
a collective dipole vibration, whereby the neutron (or 
proton) skin oscillates against an (approximately) iso- 
symmetric core. Such a collective resonance is usually 
called a pygmy resonance. Alternatively, the question 
formulated in configuration language is whether a collec¬ 
tive resonance develops, through coherent excitations of 
many valence neutrons (or protons). The answer holds 
particular relevance for validating the predictive power 
and limitations of nuclear structure models in unexplored 
regions of the nuclear chart. It is related to predictions 
about shell structure, the softness of the symmetry en¬ 
ergy, phonon coupling, and broad resonances in the con¬ 
tinuum. The reader may consult Refs. [1-3] for reviews 
on the status of the above issues in the recent past. 

The majority of studies devoted to the above is¬ 
sues, both theoretical and experimental, have focused on 
neutron-rich nuclei. One reason is that the highly accessi¬ 
ble heavy stable nuclei are, in an obvious sense, neutron- 
rich. The other side of the stability valley, comprising 
proton-rich nuclei, is no less interesting and challenging. 
Because the N = Z line reaches the proton drip line by 
Z « 50, nuclei with a few excess protons (Tz < —3) reach 
up to medium mass only. Their structure is determined 


to a large extent by the Coulomb field, besides the sym¬ 
metry energy, and of course by shell effects, pairing, and 
their proximity to the drip line. 

Early theoretical work, within the relativistic self- 
consistent quasi-particle random-phase approximation 
(RQRPA), predicted the development of a proton-skin 
oscillation in proton-rich Y = 20 isotones and very 
neutron-deficient Ar isotopes [4]. The fragmentation of 
such a resonance was further studied in Ref. [5] using 
a unitarily transformed Argonne V18 potential. Within 
the relativistic model, neutron-deficient Ar isotopes were 
predicted to develop a halo-like structure [1], character¬ 
ized by a diffuse proton distribution, favoring soft dipole 
excitations. The somewhat heavier Y = 20 isotones, on 
the other hand, were predicted to develop less diffuse pro¬ 
ton skins. Results on "^^Ni within the Skyrme-continuum 
random phase approximation predicted a similar amount 
of strength as RQRPA [6], but the origin could be single¬ 
particle excitations into the continuum. Studies of light 
Ne and F isotopes, within the RQRPA [7] and the shell 
model [8], predict an exotic proton-dominated transi¬ 
tion in ^^Ne. A Gogny-QRPA model preticts a col¬ 
lective low-lying proton excitation in ^®Ne [9]. Large- 
scale continuum quasiparticle random-phase approxima¬ 
tion (CQRPA) calculations predict up to five percentage 
points of the El energy-weighted sum to be exhausted 
by transitions below 10 MeV in various proton-rich nu¬ 
clei [10]. Recent measurements in 32,34^j, isotopes [11] 
will help assess some of the above predictions. 

In the present theoretical study we focus on Y = 20 
isotones. These nuclei are heavy enough to present 
proton-skin candidates (as opposed to proton-halo nu¬ 
clei). The Y = 20 isotones include the stable and sym¬ 
metric ^°Ca, and the magic drip-line nucleus "‘^Ni. The 
mirror nucleus of the latter, namely ^®Ca, is stable and 
thoroughly studied and can be used to pinpoint the role 
of the nucleon separation energy. 

One of our main motivations is that proton-rich nu¬ 
clei in this mass region could be accessible in future 
radioactive-ion beam facilities. Furthermore, as a testing 
ground, they can help us address a number of questions 
which have risen in past investigations. The questions 
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include, for example, the degree to which covariant and 
non-relativistic models agree in the case of proton-rich 
nuclei [6], and the role of shell effects and the contin¬ 
uum in determining the strength distribution (cf. a re¬ 
lated study of Ni isotopes [12] and the thorough inves¬ 
tigation of Ref. [13]). In the present work we use two 
non-relativistic models with complementary capabilities: 
the self-consistent quasiparticle random-phase approxi¬ 
mation (QRPA) with the Gogny DIS interaction and the 
continuum random phase approximation (CRPA) with 
Skyrme interactions. The former includes pairing effects, 
while the latter offers a proper treatment of the contin¬ 
uum. 

Our article is organized as follows. In Sec. II we intro¬ 
duce the theoretical models and the quantities of interest. 
In Sec. Ill we present our results for the N = 20 isotones. 
We summarize our results and conclusions in Sec. IV. 


II. THEORETICAL FRAMEWORK 

Here, we briefly discuss the models used in the present 
study and essential physical quantities for our analysis. 

For a study of weakly bound nuclei far from the valley 
of stability, it is indispensable to consider pairing correla¬ 
tions. A self-consistent approach of mean field and pair¬ 
ing correlations is the Hartree-Fock-Bogoliubov (HFB) 
theory which has been successfully applied not only in 
stable nuclei, but also in exotic nuclei with a large pro¬ 
ton or neutron excess. In the QRPA the input ground 
states of (open-shell) nuclei are from the HFB theory in 
the canonical basis, in which HFB ground states can be 
expressed in terms of the BCS variational parameters and 
in which the single particle density matrix is diagonal. In 
this study we work with the self-consistent quasiparticle 
random-phase approximation (QRPA) with the Gogny 
DIS interaction [14] (QRPA-I-DIS). The latter has pro¬ 
vided interesting results for the electric-dipole response 
in stable and neutron-rich nuclei [12, 15-18] which com¬ 
pare well with existing data. 

Despite its great success in describing collective exci¬ 
tations of stable and unstable nuclei, the QRPA fails to 
treat highly excited states with continuous spectra prop¬ 
erly since typical QRPA calculations are based on dis¬ 
crete spectra. This issue is especially relevant for very 
weakly bound systems such as the proton-rich ones be¬ 
ing studied here. To deal with the coupling between 
bound states and (quasi)particles in the continuum such 
as a two-particle configuration with one nucleon in the 
discrete bound level and the other in the continuum, 
one needs to construct the (Q)RPA with the appropri¬ 
ate boundary conditions and, ideally, starting from the 
ground state of a nucleus calculated using couplings to 
the continuum. A continuum RPA (CRPA) without pair¬ 
ing correlations, based on the Skyrme-Hartree-Fock cal¬ 
culations of the nuclear ground state [19] is employed in 
[6]. The same CRPA with the SLy4 Skyrme interaction 
will be employed in this work (CRPA-|-SLy4). The SLy4 


functional is chosen because it predicts similar nuclear- 
matter properties of possible relevance as the DIS func¬ 
tional, namely symmetry energy (J = 32,31 MeV for 
SLy4, DIS, respectively) and nucleon effective mass 
(to* = 0.7to) and they are both asy-soft [20]. The 
difference in the Thomas-Reiche-Kuhn (TRK) enhance¬ 
ment factor (0.25 in SLy4 but 0.6 in DIS) is expected 
to influence the peak energy of the giant dipole reso¬ 
nance (GDR) [21], but not the properties of pygmy dipole 
strength [22-24]. 

We will also discuss results with other Skyrme param- 
eterisations. 

The single-particle basis in which the HFB equations 
are solved consists of 15 harmonic-oscillator shells in this 
work. The length parameter 6 ho is such that the ground- 
state energy is minimized. The QRPA equations are 
formulated in the HFB canonical basis. By solving the 
QRPA equations, we can obtain the excitation energy Ei, 
and the two-quasiparticle amplitudes, X'' and Y‘', of the 
ly-th QRPA state. From these quantities we can calculate 
the transition probabilities and all other quantities of in¬ 
terest. For details on the QRPA implementation and its 
performance the reader may consult Ref. [25]. 

The CRPA method is formulated in coordinate space 
using Green’s functions. First the Skyrme-Hartree-Fock 
(SHF) equations are solved to obtain the single-particle 

potentials and densities. The unperturbed ph propagator 
qo 

is then evaluated, as well as the residual ph interac¬ 
tion Vi-es, in matrix form, where the rows and columns 
correspond to radial points on a mesh, (r, r'), and differ¬ 
ent blocks correspond to different operators in the resid¬ 
ual interaction and the two possible values of isospin t. 
The RPA propagator is then calculated from the 

equation 

^RPA ^ G°yres)”^G° (1) 

by matrix inversion. We use a radial mesh of 16 fm for 
and 18 fm for ^®Ni and a step size of 0.08 fm. 

In this work only spin-independent terms of the resid¬ 
ual interaction is considered in CRPA. Therefore, our 
CRPA implementation is not fully self-consistent. To 
correct for the missing effects we scale the residual in¬ 
teraction by a factor, such that the spurious state is ob¬ 
tained at zero energy. The required factors for ^*^Ca, 
48Ni^ 48ca equal 1.03, 1.11, 1.14, respectively, when the 
SLy4 functional is used. 

To investigate the properties of A = 20 isotones (and 
■^^Ca) we evaluate several quantities of interest using the 
two models. A quantity of primary interest is the transi¬ 
tion probability from the (QRPA/CRPA) excited state to 
the ground state, a measure for the dynamical response 
of a nucleus to an external operator Q. In our case this 
will be an isoscalar (IS) or isovector (IV) dipole operator. 
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The IS and IV electric dipole operators are given by 



N ■ ' Z ■ ' 

Qei = ^ rpYiMiip) “ ^ r„yiM(fn) (2) 

p—1 n—1 

in standard notation (see, e.g., Ref. [25]). In the frame¬ 
work of our self-consistent QRPA method the above E\ 
operator is equivalent to its uncorrected form, Q'p^i = 
6 X[p=i TpYiMi^p), in which only protons contribute. 

The transition matrix element of the electric dipole 
operator with spherical symmetry in the QRPA reads 

bqq' = {Xgq, - Y^g,){VqUq> - Vq>Uq){q\Q\q') , (3) 


operator 

(Ip'^(r) = (j^|^(5(r-ri)|0) . (7) 

i 

In the QRPA, the radial part of the transition density in 
the 1 “ channel can be expressed as 

6pt = {X^g, - Y^g,)(vqUg - v^uq)Rq{r)Rq, (r)(S) 

q<q',tq=t'q 

where t is for the isospin (neutron or proton) and R{r) 
is the radial wave function of the QRPA state. In the 
CRPA model, transition densities can be evaluated as 
follows. Near a resonance of interest, at energy E, the 
radial Green’s function is proportional to a product of 
radial transition densities [27], 


where u and v are the occupation probability factors of 
single-nucleon states in standard notation. The transi¬ 
tion strength B{Ey) is then given by 

2 

B{E,)= Y^q' • (4) 

q<q' 

A smooth transition-strength distribution can be ob¬ 
tained by folding with a Lorentzian, 


R{E) = Y 


■ (A - E,Y + r 2 • 


Here, 2r is the desired full width at half maximum 
(FWHM) for each isolated transition. 

We note that, in existing studies using the QRPA-I-DIS 
model [16-18], it was found possible to reproduce the 
measured El strength in various nuclei by considering a 
3 MeV downward shift, i.e., by accepting a global dis¬ 
crepancy of 3 MeV in the low-energy region between 
the model and the data. The origin of this systematic 
discrepancy is not clear, but it may lie in convergence 
issues or the need to include effects beyond RPA (cf. 
Refs. [13, 26]). In the discussion of our results it is made 
explicit when we take into account this empirical shift. 

In the CRPA the strength distribution is determined 
by the polarization diagram via the Green’s function 


R{E) = = ““51 / d^r<fr'ft{r)Gtt'{r,E]E)fti{r), 

4 . 4 ./ 'J 


where /t(r) denotes the spatial and isospin dependence 
of the external field. The obtained R{E) automati¬ 
cally includes the escape width. The distribution can 
be smoothed further for presentation purposes by adding 
an imaginary part to the ph energy in G°, analogous to 
the parameter T in the Lorentzian of Eq. (5). 


The transition density is also an essential quantity to 
study the dynamics of nuclear collective motion, and it 
is defined as the transition matrix element of the density 


9Gt,t/(r,r'; A) oc 5pt{r-,E)5pt'{r'\E). (9) 

Folding the above quantity with the radial part of the 
external field, /t(r), we get a function 

5p[^\r-E) =Y f ^GtAr,r'-,E)fAr')dr', (10) 

t' d 

which is proportional to 6pt{r; E), if Eq. (9) holds. If we 
now fold this with ftir), we obtain tt times the transi¬ 
tion strength per energy unit. Although the above fac¬ 
torization is not always perfect, the procedure helps us 
to obtain a quantity, loosely called a transition density, 
which determines how the nuclear density couples to the 
external field at a given energy. 

The last quantity that we will examine in this work 
is the energy weighted-sum (EWS) of El strength. In 
QRPA we may write it as 

mi(A,)= E,-B{Ei), (11) 

i:Ei<Ej 

By m\°^ we will denote the sum over all the available 
QRPA states. In CRPA the equivalent quantities are ob¬ 
tained by an energy-weighted integration of the transition 
strength probability. 


III. RESULTS 

From the shell-structure point of view, the N = 20 iso¬ 
tones differ only in the occupation number of the proton 
0 / 7/2 orbital, which is zero for ^°Ca and 8 for "‘^Ni. The 
orbital can be spatially extended if its binding energy is 
small. According to the HFB model with the DIS func¬ 
tional, that energy is lower than 1 MeV in ^®Fe and ^®Ni. 
According to the SHF model with the SLy4 functional, 
the orbital is somewhat unbound in ^®Ni, but localized 
via the Coulomb barrier. This is in contrast to the well- 
bound neutron O/7/2 orbital in all mirror nuclei. Ac¬ 
cordingly, the proton separation energies in the N = 20 
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FIG. 1. (Color online) Radial wavefunction of the proton 
0 / 7/2 orbital in ^®Ni and of the neutron O/ 7/2 orbital in '^^Ca 
according to the Hartree-Fock model with the SLy4 functional 
(solved in coordinate space) and the DIS functional (in a 
harmonic-oscillator basis). The inset shows the same in semi- 
logarithmic scale. 

isotones are quite different from the neutron separation 
energies in the mirror Ca isotopes. From the binding en¬ 
ergies tabulated in the 2012 evaluation [28] we find that 
the neutron separation energies in the Ca isotopes well 
exceed 9MeV. For the Z = 20, 22,..., 28 isotones, on 
the other hand, the proton separation energies take the 
values 8.32, 3.75, 3.03, 1.54, 0.692MeV, respectively. 

In Fig. 1 we show the corresponding radial wavefunc- 
tions in '^®Ni as well as the neutron O/ 7/2 wavefunction of 
the mirror nucleus ^®Ca according to both models. Both 
models produce a more-extended proton wavefunction in 
^®Ni than the (much more bound) mirror neutron wave- 
function. This result is dictated by basic quantum me¬ 
chanics as it originates in the much weaker binding of the 
orbital in the proton-rich counterpart than the neutron- 
rich one. Also visible is the different asymptotic behav¬ 
ior predicted by the two models. The extended tail of 
the wavefunctions produced by the Hartree-Fock model 
solved in coordinate space is absent form the results of 
the model solved in the harmonic-oscillator basis. 

The difference between the proton and neutron root 
mean square radius, — Rn, defines the thickness of 
the proton skin. The absolute values obtained with the 
HFB-I-DIS model are listed in Table I. Also listed are 
the absolute values for namely the mirror nu¬ 

clei of ^®Fe and ^®Ni, corresponding to the neutron-skin 
thickness. The values |i?p — i?„| for '*°Ca, “^^Ca and '^^Ni 
obtained with the SHF-|-SLy4 model are, respectively, 
0.076, 0.111 and 0.281 fm. We observe that the pro¬ 
ton skins of the proton-rich nuclei are much thicker than 
the neutron skins of the mirror nuclei. Therefore, even 
though no skin-oscillation has been observed in "^®Ca at 
low energies [18], it remains possible that a proton-pygmy 
mode develops in "‘^Ni or '*^®Fe - a possibility which we 
aim to assess in this work. 

In Fig. 2 we compare the GDRs obtained from the 
two models and from experiments [29], if available. The 
CRPA model can be applied for closed-shell configura- 
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FIG. 2. (Golor online) Isovector dipole strength distributions 
smoothed with a Lorentzian of F = 1 and experimental data 
where available. 

tions only, therefore no results for ^®Fe can be shown. 
The discrete isovector dipole strength distributions are 
smoothed by a Lorentzian of width F = 1. Roughly 
speaking, the peak position of the distributions from the 
two models differ by a few MeVs, which can be attributed 
to the different TRK enhancement factors. For ^°’'^®Ca, 
it is seen that the observed GDR peak lies between the 
predictions of the two models; CRPA-|-SLy4 describes the 
low energy region of the GDR fairly, while QRPA-I-DIS 
accounts for the tail of the GDR at high energies better. 
In this sense both models are found equally appropriate 
in the case of the GDR. The fragmentation patterns are 
similar with both functionals. 

We stress that, athough the theoretical predictions for 
the GDR peak have been found correlated with the pre¬ 
dicted enhancement factor, no correlation has been found 
between the pygmy dipole strength and the enhancement 
factor [22-24]. 

Next, we examine in detail the low-energy transitions 
with model QRPA-t-DlS and then we compare with 
CRPA results for further insights. 


A. QRPA results 

Figure 3 shows the IS and IV dipole transition strength 
from QRPA calculations for the N = 20 isotones and 
for ^®Ca. In the IV response the main feature in all 
nuclei is the GDR. In the IS response we observe at 
least two dominant features: The concentration of high- 
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lying peaks, corresponding to the dipole compression 
mode [30], and the strong transition around 10 MeV 
(9.9 MeV < E < 10.8 MeV). The latter corresponds 
to the IS surface dipole mode, or IS low-energy dipole 
(IS-LED) mode, studied also in Refs. [12, 15-17, 31]. 
All structural properties of the IS-LED calculated in the 
present work are found in line with those studied in the 
above-cited publications, to which we refer the reader for 
details. As an example, we show in Fig. 4 the transition 
densities of the IS-LED peaks in "^^Fe and ^®Ni. The pro¬ 
ton and neutron transition densities are of the IS-LED 
type, namely in phase and with a node near the surface. 
We have checked that the corresponding transition den¬ 
sities in the other nuclei examined here exhibit a similar 
trend to the ones in Fig. 4. In addition, the IS-LED 
peaks account for very small E\ strength as manifest in 
Fig. 3 and Table 1. The corresponding states have been 
observed in around 7 MeV [18, 32] and other nu¬ 

clei [31]. We point out that in ^®Fe, "^^Ni the IS-LED 
would be unbound because of the low particle emission 
threshold. 

In Fig. 3 we observe also the effect of adding more 
and more protons to the symmetric nucleus ^*^Ca, namely 
increased transition strength, in both the IS and the IV 
channel, at energies between the IS-LED and the GDR. 
Nuclear properties that change as we add protons, are 
the asymmetry |A^— Z\lA, the occupation number of the 
valence ttO/y /2 orbital, and the proton separation energy. 
The effect of the latter can be isolated by a comparison 
of ^®Ca and "‘^Ni. Indeed, if we compare the spectra 
of these two mirror nuclei and the values of low-energy 
strength tabulated in Table I, we find that the barely 
bound ^®Ni is much more strongly excited than the stable 
^®Ca. Similarly, if we compare the values of low-energy 
strength tabulated for the A = 46 mirror nuclei, we find 
that "^^Fe is much more strongly excited than the stable 
^®Ca at low energies. In the stable Ca isotopes, the E\ 
strength summed up to 15 MeV is comparable to the 
strength of a single-particle transition. 

The fragmented resonant-like structure we have ob¬ 
tained between the IS-LED and the GDR in the proton- 
rich isotones represents a candidate for a proton-skin os¬ 
cillation. Next, we assess this possibility. First we note 
that a skin oscillation would respond strongly in both 
channels, IS and IV. This is in analogy to the strong 
IS response expected of neutron-skin modes [33-35]. In 
Fig. 5 we display the transition densities of states which 
appear relatively strong in both channels, for the proton- 
rich isotones ^"^Gr, ^®Fe, ^®Ni. We have verified that the 
transition densities of omitted states do not show any 
skin-oscillation pattern or they resemble those of neigh¬ 
boring states shown in Fig. 5. The states beyond 15 MeV 
(rightmost states in Fig. 5) show a strong IV character 
and can be considered as fragments of the GDR. Gener¬ 
ally, the leftmost state for each isotone appears the most 
promising proton-pygmy candidate. In order to assess 
its collectivity we examine which quasiparticle configura¬ 
tions contribute to the El strength. In Fig. 6 we show the 


quantity 5^^,, Eq. (3), for the IV operator Qei- We notice 
many contributions of destructive coherence in ^®Fe, even 
if we consider only proton configurations, i.e., the un¬ 
corrected operator (see below Eq. (2)). In neutron- 
rich nuclei, neutron-skin oscillations have been associated 
with destructive coherence [34, 36]. In ^®Ni the dominant 
configuration is 7r0/7/2 —> 'Kld ^/2 (tt stands for protons). 
In "^^Gr there is also a strong 7r0/7/2 —>■ 7 r 059/2 compo¬ 
nent. In "^^Fe, however, the transition is evidently not of 
single-particle character. The value in the third column 
of Table I in comparison to the first column support the 
conclusion that the mode in question is stronger than a 
single-particle transition. 

We can conclude, based on our QRPA-fDlS results, 
that we expect a proton pygmy mode to develop in 
proton-rich = 20 isotones, most unambiguously in 
"‘^Fe. Based on Figs. 3, 5 we expect a proton-pygmy 
resonance to appear somewhat fragmented and possibly 
surrounded by other El excitations. 

We note that the relativistic models of Refs. [4, 6] 
predict, as lowest-lying states, proton-pygmy modes of 
enhanced El strength. By contrast, the QRPA-I-DIS 
model predicts an ordinary IS-LED of very little El 
strength. The situation is reminiscent of that in neutron- 
rich nuclei, where neutron-skin vibrations are predicted 
for rather moderate neutron excess [12, 17]. In proton- 
rich nuclei the two kinds of states may coexist in the 
continuum in reality. Next we address the treatment of 
the continuum. 


B. Comparisons with CRPA results 

So far we have analyzed the QRPA eigenstates as di- 
crete and localized excitations. There are two effects 
missing from such a treatment, related to the fact that 
the particle-emission threshold is very low in the proton- 
rich isotones. First, practically all excited states are un¬ 
bound and they may have a finite escape width. Second, 
protons can be excited to non-localized states near the 
Goulomb barrier (and neutrons to states of small positive 
energies), which cannot be efficiently described in a har¬ 
monic oscillator basis. The effect of the basis was already 
evident in the Hartree-Fock results of Fig. 1. The CRPA 
can address the above-mentioned effects, albeit only for 
closed-shell configurations. Next we examine our CRPA 
results on ^®Ni and its mirror nucleus ^®Ca in comparison 
to our QRPA results. 

We begin by comparing the response of the two nu¬ 
clei, in both models, in Fig. 7 and in Tables 1,11. At 
this point we are discussing the top two rows of Table II 
(SLy4 functional). Each model predicts similar patterns 
for the GDR of the two nuclei, but enhanced low-energy 
El strength in ^®Ni compared to ^®Ca. In the case of 
CRPA, we recognize that the distribution shows a large 
escape width, as reported also in Ref. [6]. 

In Fig. 8 we compare the results of the two models 
for each nucleus in the low-energy regime. There is a 
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FIG. 3. (Color online) QRPA IS and El dipole strength distributions in units of e^/m® and e^/m^, respectively. Note that 
the IS strength is divided by 250 for presentation purposes only. 


TABLE I. Single-particle unit (single-proton for N = 20 isotones and single neutron for ^®’‘^*Ca), El transition strength of the 
IS-LED, and El strengths summed up to (13, 15, 17) MeV, in units e^fm^ calculated with the QRPA-I-DIS model. The proton 
(or neutron) skin thickness is also given in fm. 



Bs.p 

B{Ei) 

E B{E^) 

Ei<13 

E b{e,) 

Ei<15 

E B{E,)\Rp-R^\ 

E,<17 

40 Ca 

0.188 

0.004 

0.008 

0.117 

0.179 

0.042 

42Ti 

0.177 

0.008 

0.116 

0.249 

0.361 

0.106 

44cr 

0.166 

0.009 

0.189 

0.320 

0.576 

0.159 

46Fe 

0.156 

0.008 

0.251 

0.421 

0.835 

0.207 

48Ni 

0.148 

0.009 

0.161 

0.553 

1.211 

0.249 

46Ca 

0.156 

0.003 

0.050 

0.149 

0.443 

0.108 

48Ca 

0.148 

0.006 

0.006 

0.089 

0.763 

0.145 


large amount of El strength predicted by CRPA in 4®Ni 
below 10 MeV, corresponding to more than 1% of the 
energy-weighted sum, but not in As evident from 

the values in Tables I, II, this amount of strength is not 
described by QRPA-I-DIS, unless we take into account 
an energetic shift larger than the 3 MeV introduced in 
Sec. II. 

In Fig. 8 we include the results of the CRPA-|-SLy4 
model with box boundary conditions imposed, to elimi¬ 
nate the escape width. For 48Ca the effect is very weak. 


For 48Ni the strongest effect is in the narrowing of the 
low-energy tail into three peaks, indicating that there are 
at least three configurations contributing to this broad 
structure, but the total strength does not change. 

The two models, namely DIS and SLy4, have similar 
symmetry properties (the SLy4 functional is even more 
asy-soft than the DIS) and the same nucleon effective 
mass. One might expect more strength predicted by 
the Gogny DIS functional, but we observe the opposite. 
Therefore, the only factor we can identify as enhancing 
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46 Fe *»Ni 



FIG. 4. (Color online) Proton and neutron transition den¬ 
sity of isoscalar low energy dipole states of ‘*®Fe and ^®Ni at 
E=10.48 MeV and E=10.35 MeV, respectively. 
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FIG. 5. (Color online) Proton and neutron transition densities 
for selected isotones. 


the low-energy strength in the case of CRPA-|-SLy4 is the 
proper treatment of spatially extended configurations. 
The QRPA cannot account for the low-energy transi¬ 
tion strength because of the bad convergence properties 
of the harmonic-oscillator basis for extended wavefunc- 
tions. Our statement is corroborated by a comparison 
of the QRPA results using different length parameters 
6 ho- a. larger value of 6 ho (1-90 fm instead of 1.65 fm 
in '‘^Ni and 1.60 fm in ^®Ca) has the effect of shifting 
pygmy strength to lower energy (and a rather weak over¬ 
all effect at high energies) in “^^Ni. The effect is weak 
but present in '^^Ca. The sensitivity of the results to the 
length parameter signifies that a proper treatment of the 



E„. [MeV] [MeV] [MeV] 


FIG. 6. (Color online) Two-quasiparticle configuration contri¬ 
butions to the transition matrix element of the pygmy states 
of '‘■‘Cr, ^®Fe, and '‘®Ni. 


TABLE II. El transition strength strength summed up to 
(10, 13, 15) MeV, in units e'^fm^, and percentage points of 
the energy weighted sum below 10 MeV. 



/ dBiE) 

E<10 

/ dB{E) 

E<13 

lOOm'T"" 

m\°^ 

SLy4:48Ni 

0.242 

1.241 

1.01 

48 Ca 

0.028 

0.560 

0.11 

BSk5:48Ni 

0.254 

0.950 

1.08 

48 Ca 

0.029 

0.528 

0.12 

SkI3:48Ni 

0.517 

1.277 

2.11 

48 Ca 

0.080 

0.564 

0.37 


threshold strength is very important in the description of 
pygmy proton strength. 

Finally, we discuss Fig. 9, which displays the cummu- 
lative energy-weighted sum of El strength in both mod¬ 
els, the total being normalized to one. According to the 
CCRPA+SLy4 model, 1% of the energy-weighted sum 
can be expected below 10 MeV in ^®Ni, in stark contrast 
to the mirror nucleus and to ^°Ca. The importance of the 
low separation energy, rather than just the asymmetry, in 
generating this amount of strength is evident. According 
to the QRPA+DIS model, as long as we take into ac¬ 
count an energetic shift of about 3 MeV, for the reasons 
discussed in Sec. II, a similar percentage of the m\°^ is 
predicted in all isotones with Z > 20, with the expection 
of ^®Fe, where this amount may even double. (Without 
the 3 MeV shift the predicted amount of strength evi¬ 
dently vanishes.) As we conclude from the values in Ta¬ 
ble I, 1% is consistent with one single-particle unit. '^^Fe 
is therefore the most promising candidate for detecting 
a coherent skin mode. Nonetheless, a measurement in 
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FIG. 7. (Color online) The El strength distributions for and ''^Ni in the two models. A Lorentzian with F = 0.2 is 
employed for the QRPA+DIS results. 


more-accessible lighter isotones would also be of value: 
If stronger transitions are detected, they might signify 
a collective effect. We note that our prediction for 1% 
of the EWS should be considered as a lower theoretical 
limit since it comes from asy-soft functionals with a low 
effective mass. 

In the above we have focused on the isovector, El tran¬ 
sition strenght. As is evident from Fig. 3 and our results 
with the Skyrme functional. Fig. 8, the pygmy resonance 
and other proton transitions would be very strongly ex¬ 
cited by isoscalar probes. An issue to keep in mind in that 
case is that a proton-skin mode and the ordinary IS-LFD 
mode would both be unbound in proton-rich nuclei and 
likely indistinguishable. Helpful in this respect should 
be that the IS strength of the IS-LFD transition in all 
isotones and in the mirror stable Ca isotopes should be 
approximately the same. Therefore, it can be considered 
approximately known, because the strength has been ex¬ 
tracted in 32] and compares well with our 

calculations. A large excess of IS strength in proton-rich 
nuclei could then be attributed to exotic effects. 


C. Other functionals 

We have examined CRPA results with two more func¬ 
tionals and the results are tabulated in Table II. The 
BSk5 [37] functional has similar symmetry-energy prop¬ 
erties as the DIS and SLy4 functionals [38], but a higher 
nucleon effective mass, namely 0.92 m. The isovector ef¬ 
fective mass is also somewhat higher for BSk5 (0.84m) 
compared to SLy4 (0.80m). The GDR is found some¬ 


what shifted to lower energies in comparison to the SLy4 
prediction, but almost the same El strength distribution 
is predicted below 10 MeV. 

In Ref. [6] it was observed that a Skyrme model of 
high isovector effective mass m* is required to obtain 
more strength and thus agreement with the relativistic 
models used in that work. The functional MSk7 with 
m* = 1.05m predicted the largest amount of low-energy 
strength, even though it is among the most asy-soft used 
there. We now test whether a combination of a low m* 
and a high L can also give a large amount of strength. To 
this end we use SkI3, which predicts a low isoscalar effec¬ 
tive mass but L = 100.53 (the DD-MEl and DD-ME2 rel¬ 
ativistic models of Ref. [6] predict L = 55.42, 51.24 MeV, 
respectively [20]) and an isovector effective mass simi¬ 
lar to the other models, namely 0.8m. The results are 
given in Table II. The amount of strength below 10 MeV 
is larger than twice the amount predicted by the other 
functionals. We therefore postulate, based on this ex¬ 
ploratory calculation, that the asy-softness of the func¬ 
tional can affect the amount of pygmy dipole strength 
in an analogous way as the neutron-pygmy strength (see 
also Ref. [39]) and that both the isovector effective mass 
and the slope parameter can enhance the pygmy-dipole 
strength in proton rich nuclei. In Ref. [6] mainly the 
former is presumably at work in the case of the MSk7 
functional, and the latter in the case of the relativistic 
models. Obviously, further systematic studies must be 
undertaken, along the lines of Refs. [23, 40-42], before 
the above trends can be established. 

In beyond-mean-field models, an enhancement of the 
effective mass can be implicit through coupling of single- 
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FIG. 8. (Color online) The low-energy El (or IV) and IS 
strength distributions for and ‘*®Ni in QRPA-fDlS, 

QRPA-fDlS with a larger length parameter 6 hO) Hie 
CRPA-|-SLy4, and the CRPA-|-SLy4 with box boundary con¬ 
ditions. A Lorentzian with P = 0.2 is employed. 




E[MeV] E[MeV] 


FIG. 9. (Color online) The energy-weighted sums of Af = 20 
isotones and ^®Ca in the two approaches. 


particle states to surface vibrations. 


IV. SUMMARY AND CONCLUSION 

We have studied theoretically the electric dipole re¬ 
sponse of the V = 20 isotones at low energies (mainly 
below 10 MeV). We analyzed results of the QRPA solved 
in the harmonic oscillator basis and using the Gogny 
DIS interaction and of the CRPA solved in coordinate 
space using Skyrme functionals. These are widely used 
phenomenological methods, which in the future shall be 
complemented by more ab initio approaches [43]. 

Larger amounts of El strength in the asymmetric 
N = 20 isotones are predicted than the amounts of 
strength predicted or detected in equally asymmetric 
Z = 20 mirror nuclei. The difference between the two 
mirror sets of asymmetric nuclei, Z = 20 or N = 20, is 
the much looser binding of the latter. Thus our results 
point unambiguously to the important role of structural 
effects, as opposed to global parameters like asymmetry, 
in determining the El spectrum at low energy. 

An exotic collective excitation is found most likely in 
"‘^Fe but perhaps also in "^^Cr and '*®Ni, exhausting more 
than 1% of the energy-weighted summed strength. A 
similar amount of strength was predicted in Ref. [10] 
for the N = 20 isotones using the relativistic contin¬ 
uum QRPA with a relatively soft relativistic model. In 
our above prediction we take into account the results 
of the CRPA solved in coordinate space, the expecta¬ 
tion that the energies of the low-lying states are over¬ 
estimated by the QRPA model by about 3 MeV [12, 16] 
and that the treatment of extended wavefunctions by the 
harmonic-oscillator basis in the same model is inefficient. 
A correct, converged treatment of threshold transitions 
(as in CRPA), and therefore of extended wavefunctions, 
is found important for the description of proton pygmy 
states since the proton emission treshold is extremely low. 

We consider the predicted amount of low-energy 
strength, namely 1% of the EWS, as a lower theoreti¬ 
cal limit, because it comes from asy-soft functionals with 
a low effective mass. We predict furthermore that the 
strength is likely distributed in a broad structure in the 
spectrum, which would respond very strongly to isoscalar 
probes. 

Based on an exploratory comparison of different func¬ 
tionals, we postulate, pending systematic studies, that 
the effective mass as well as the slope parameter L can en¬ 
hance the low-energy strength. A systematic and quan¬ 
titative investigation of the two effects for proton-rich 
nuclei, along the lines of similar studies in neutron-rich 
nuclei [23, 40-42] and with proper convergence, is highly 
desirable. 
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